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Abstract 

Incommensurate easy flow channels in an otherwise perfect vortex lattice are 

investigated. The associated (point) defects in the lattice inside the channel 

cause an almost vanishing critical current, as shown by molecular dynamics 

simulations and a comparison with the Frenkel-Kontorova model. In addition 

to the normal flux flow behavior we find a low mobility regime at small drives 

associated with defect motion. We treat this situation analytically for the 

case of a single defective vortex row. We also briefly discuss the relation to 

existing experiments on artificial vortex channels. 
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In recent years, vast interest has grown in the dependence of static and dynamic properties 
of vortex matter on the spatial profile of the pinning potential originating from defects in 
the underlying material. The behavior of the vortex lattice was investigated experimentally 
as well as theoretically in the presence of a wide variety of these pinning landscapes, rang- 
ing from purely random to highly periodic [Q. Several of these studies revealed the 
existence of static channels of easy vortex flow at currents just above the critical value [@|§. 
In addition, this critical current was shown to be proportional to the interaction strength 
between vortices inside the channel and those in the channel edges (CE's), as expressed by 
the shear modulus of the vortex lattice cqq A precise theoretical description of this 

phenomenon and of its dependence on the channel width is still lacking. In this paper we 
study for the first time the properties of static easy flow channels within the framework of 
the Frenkel-Kontorova (FK) model ||. It is demonstrated that a mismatch between chan- 
nel width and lattice constants induces (point)defects in the channel leading to an almost 
vanishing shear strength. This can have important implications for the properties of vortex 
matter in a pinning potential with large spatial variations in strength. In addition, our re- 
sults on the dynamics of incommensurate structures may be applied to various other fields, 
such as vortex dynamics in overdamped Josephson junction arrays, transport properties of 
charge density waves [0 and solid friction of confined layers |§||. 

We consider 2D vortices at T = in an easy flow channel of width w and length /, 
confined by two semi infinite pieces of a rigidly pinned vortex lattice with perfect hexagonal 
structure of lattice constant a and row spacing b = a \/3/2. The pinned structures forming 
the edges are chosen to have their principal axis in the x-direction parallel to the channel 
and have a relative shift Ax (left upper corner of Fig. la). First we describe qualitatively 
the case in which w ~ b and Ax = (restricting the degrees of freedom for the mobile 
vortices to the x direction), corresponding to the single chain FK model. Molecular dynamics 
(MD) simulations will then provide a justification of these qualitative arguments as well as 
show that our line of reasoning also applies to larger channel widths and a relative CE-shift 
Ax + 0. 
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When a shear force is applied to the vortices in the channel, the onset of flux flow 
appears as the force exceeds the maximum shear stress. The maximum shear force density 
F s may be obtained from c§§: F s = 2Ac§§jw where A accounts for the lattice anharmonicity 
and orientation of the shear direction with respect to the principal lattice axes |T0| . For 
our orientation and a harmonic potential from the vortices in the CE's, it follows that 
A = Then, for w ~ bo, the 'edge' potential (per unit length) for a single vortex 

inside the channel can be defined as U — |(1 + cos^o^)), with a = a 3 C66/(vr 2 6 ) and 
ko = 2n/aQ. Next, we consider the mutual interaction between the mobile vortices in the 
channel. These vortices have an average spacing a which is given by the induction B as 
a = §o/(Bw) = aobo/w with $ the fluxquantum. We use a harmonic approximation for the 
interaction, while its nonlocal nature is preserved by including interactions with n neighbors 
in the chain. The effective spring constants K n then can be obtained from a London 1/r 
interaction force between vortices: K n = (§l/4:7tfio\ 2 )(l/na) 2 with A the penetration depth. 
Thus, in the presence of a transport current density J applied perpendicular to the channel 
, the equation of motion for vortex i in the channel reduces to an overdamped FK model: 

lii = f + iism(k Xi) + K n (x i+n + Xi- n - 2xj), (1) 

n 

where / = J$o is the Lorentz force, /i = ira/ao is the amplitude of the periodic force, and 
the friction coefficient 7 is related to the flux flow resistance by 7 = BQ / pjf(B). 

In the FK-model the critical force f c and the dynamic properties largely depend on both 
the value of a/a as well as the dimensionless elastic constant g = a 2 ) J2n n2 Kn/2'K 2 c ! ' |7|,|TT 



From the above given expressions for a and K n and using the low field expression for the 
shear modulus cqq = $o-B/(167r/i A 2 ), we find that g ~ (A/a)(ao/a) 2 for our system. Since in 
practice [| A/a > 1 we restrict ourselves to the case g > 1 |T2| . In case a/a = 1, the vortex 
row locks in with the periodic potential and f c = fi. For rational a/ao, the row contains 
point defects with a size Id ~ 2^a^fg. Then f c is determined by their Peierls-Nabarro (PN) 
barrier |TT]] which, for an isolated defect and g > 1 is given by f c ~ (S/3)7i 3 gfie~ 7T2 ^ . 
Accordingly, f c /n 1 and the channel has an almost vanishing critical current. Finally, for 



irrational a/ao and n = 1, Aubry has shown that a transition exists to a state with f c = as 
g exceeds a critical value g c . For a/a = (V5 — l)/2, <7 C = 1 [0- The condition g > 1 then 
implies that the vortex channel can have a truly vanishing critical current. The influence of 
nonlocal elastic constants WM on the above results can be neglected as long as Id > A, i.e. 



for A ^ 50ao- 

In order to check the above scenario, we performed MD simulations of easy flow channels 
in a thin film. The dynamics of a vortex inside the channel is governed by the following 
equation: 

r^J + ^f^), (2) 

with a as unit of length, $Q/(27r/XoA 2 ao) as unit of force (per unit vortex length) and 
47r/i A 2 /(\/3pj-/) as unit of time. The external drive / is directed along the channel. The 
interaction force is chosen as f v = (l/r)(l — r 2 /r 2 ) 2 with r c = 3.33 IJ. This interaction 
results in a shear modulus coinciding with the above given expression (in dimensionless units 
c 66 = 1/(4^/3)). In Eq. (fj), the index j runs over all vortices, including those in the CE's 
(see Fig. la). We adopted periodic boundary conditions with I large enough, such that our 
results did not depend on this length scale. For each w, we relaxed the system to its ground 
state (GS) and then measured the force-velocity (f-v) characteristic by stepwise cycling the 
force (/ = — > f max — > 0) and taking the stationary velocity v = (x)i ; t at each step. 

Fig. la shows the resulting f-v curves for Ax = and w ~ bo. For w = bo, the maximum 
shear strength is observed with a value f c = 0.054, coinciding with the dimensionless value 
of /Lt in (|TJ) . Additionally, the dynamics are identical to that of a single particle in a sinusoidal 
potential. Both observations are in agreement with the fact that the 'edge' potential resulting 
from our vortex interaction is indeed harmonic. More interesting is the behavior in the 
incommensurate case (w ^ b ). As seen in the inset f c j\i ^ 0.001, which can be understood 
from our previous arguments when we calculate g: expanding /„ in (|2|) we obtain K n = 
(na/a Q )~ 2 + 0.18 — 0.003(na/a ) 2 , resulting in g 6. A typical example of the GS defect 
structure is shown in Fig. lb, where the open symbols represent the vortex density along 



the channel for w = 0.97&o- 

Turning to the dynamics in Fig. la, it is seen that the defective vortex rows posses a 
regime of low mobility in the f-v curves. In this regime, transport in the row is carried by 
the defects which have a velocity Vd much larger than the average velocity v = CdVd, with 
Cd = |1 — (ao/a) | the defect density. A snapshot of the moving row then resembles the GS 
defect structure (see again the open symbols in Fig. lb). A useful quantity here is the 
mobility of an isolated defect. From the simulation we get Vd/f ~ 13, a value which one can 
also obtain from the analytical treatment described below. The mobility of an interstitial 
is slightly larger than that of a vacancy due to the anharmonic vortex interactions. When 
entering the high mobility regime where / > //, the defects are smeared out in a slight 
sinusoidal density modulation and the transport becomes coherent within the row (see the 
filled symbols in Fig. lb). 

It turns out that we can capture the transport characteristics within a perturbation 
theory as used in Ref. Consider Eq. ([[]). It is known that the motion of all particles 
in the FK chain can be completely described by the dynamic hull function h(ia + vt) = 
Xi(t) — (ia + vt) = h((p), which represents the deviations of the particle positions from the 
undisturbed sliding values. From (]!]) we obtain the following equation for h(<p): 

Vu(l+h'(<p)) =f + fi sin(k (p + h(ip))) 

+ J2K n (h(tp + na) + h((p-na)-2h(tp)). (3) 

n 

When mapping (p back onto one period of the potential, <p = (ip — ido), h(<p) has the 
periodicity a resulting in: 

f = TV + — r\h'm 2 d(p. (4) 
a Jo 

The shape and amplitude of h(0) represent the defect structure of the moving FK chain 
and thus depend on a/ do, g and v. For large speed or in case of both large lattice mismatch 
and g ^> 1 (above Aubry's transition), h((p) is sinusoidal with a small amplitude ||. In 
our case, h{(p) can have a more complex shape as well as a large amplitude. We introduce 
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higher orders in the perturbation to describe the dynamics for any a/ao by solving (0) with 
a trial hull function of the form h{(p) = J2m=i [Hm^ tkorn ^ + c.c.J • Inserting this hull function 
in (||), we obtain a set of coupled equations for the coefficients H m . Solving for H m yields 
the approximate form of h(<p). A full account of this work will be published elsewhere. Now 
the expression describing the f-v characteristics is obtained by using this hull function in 
Eq. (H), which results for M = 3 in: 

27rV 7 t> 

where aoVL 2 c represents the elastic restoring force on nonlinear excitations in the FK chain. 
From the perturbation approach Q c is expressed in /i and the restoring forces on linear modes 
(phonons) at wavenumber k = mk : 



n c (a = a b / W ) = mWo + fiH2 > (6) 



in which fij^ = Q 2 (mk ) = ^ n 2K n (l — cos(nmafco)). The last term in Eq. (||) corresponds 
to the sliding friction in the field of tribology ||. It results from the dissipation of the 
internal (nonlinear) collective modes in the chain. The crossover from low to high mobility 
occurs at v c = aoQ 2 ./ (27ij) where the amplitude of the hull function drastically decreases. 
In order to compare (|5]) with the simulated data, we use the previously obtained expression 
for K n and take into account that fi slightly depends on w as well. The corresponding f-v 
curves, shown in Fig. la, agree well with the simulated data. The transport of a defective 
vortex row can thus be accurately described by perturbing the harmonic FK model. 

In order to study channels of larger width, we adopted an edge shift Ax(w) with a 
saw tooth shape (0 < Ax < ao/2) which assures that, as we vary w, a perfect hexagonal 
structure is retained for w = pbo with p an integer. However, for w ^ pbo, the qualitative 
behavior did not depend on Ax. The inset to Fig. 2 shows the simulated f-v curves for 
3.5&o < w < 5&o- The features are similar to the f-v curves for w ~ bo and these are again 
associated with the presence of defects in the channel. Two characteristic structures can be 
identified: for | (w/b ) — p |< 0.4 (I) the structure consists of p vortex rows plus interstitials 
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(w/b > p) or vacancies (w/b < p). The defects are not necessarily equally distributed 
within a row. In addition the defect density differs per row due to the repulsive interaction 
with the channel edges. The latter is illustrated in Fig. 3a, where two snapshots of the 
motion for w = 4.08&o in the low mobility regime are displayed with interstitials indicated 
by arrows. For | (w/b ) — (p + 1/2) |< 0.1 (II), clustering of vacancies (w/b > p + 1/2) or 
interstitials (w/b < p+ 1/2) occurs. As a result, alternating regimes form along the channel 
of p and p + 1 rows, which we refer to as stacking faults (SF's). The regimes are separated 
by dislocations with Burgers vectors at angles of 60° with the channel direction, as shown 
in Fig. 3b. 

Structure I is just a 2D extension of our results for a single incommensurate vortex row. 
The shear strength is almost vanishing due to the small PN barriers of the point defects in 
each row. The low mobility regimes in the f-v curves once more originate from motion of these 
defects. However, two additional features appear. First, recalling our simple description of 
defect motion in a row {y r = c^rVd), we see that, for equal defect velocities in different rows, 
a different defect density causes different average row velocities v r , i.e. plastic motion within 
the channel. This clearly shows up in Fig. 3a, where the interstitials move at the same 
velocity Vd v r while the row velocities vary from v 2 — to v 3 — 2v. Second, occasionally 
a hysteresis is observed in the transport curves of I (see inset to Fig. 2). This hysteresis 
originates from either a redistribution of defects between rows (for 0.2 < \w/b — p\ < 0.4) 
or within the rows. After the transition, defects within the row have the maximum possible 
separation and they are clustered with those in neighboring rows, as expected from the 
standard FK model. Now the mobility of a defect moving along a locked (defect free) region 
of a neighboring row is small compared to that of a defect traveling along with a neighboring 
defect, since part of the periodic potential is destroyed in the latter case. This explains the 
transport hysteresis and the fact that the mobility is always larger after the transition. 
However, the nature of the transition and whether or not it occurs depend both on Ax and 
w. In the high mobility regime of I, vortex transport is no longer carried by the defects 
but by the rows as a whole. They attain the same velocity and the flow loses its plastic 
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character. 

The SF's of II can locally form a commensurate structure. Then one naively expects f c 
to be of the same order as the commensurate value. The actual value, f c ~ fibo/(15w), is 
much lower because the SF's can move plastically by continuously creating and annihilating 
the dislocations at their borders. This is supported by the neighboring, defective rows: these 
rows 'feed' the SF with point defects so that the longitudinal vortex displacements change 
into transverse displacements when entering the SF, as shown by the trajectories in Fig 3b. 
Now the naive estimate, f c ~ fibo/(2w), marks the crossover to the regime where the SF's 
move coherently with the rest of the structure: vsf = v. The hysteresis, observed in the f-v 
curve for w = 4.546 , is now caused by transitions in which the number of SF's decreases. 
For / I, this appears as a reduction in the amount of plastically moving SF's and thus an 
increase in mobility. 

The main plot of Fig. 2 summarizes the behavior of f c vs. w, where we took the 
maximum f c in case of hysteresis. The continuum result f c = fib /w applies only for w = pb 
while f c of the incommensurate channels reflects our above discussion of the barriers of the 
corresponding defects. To make a link to the experiments on samples with artificial vortex 
channels as used in Ref. |J, we also show the measured shear strength f c = 27iJ c fi \ 2 a /§Q 
versus w/bo of such a device. It is seen that the commensurability peaks are lower and 
smeared out. Preliminary results of simulations in which positional disorder of vortices in 
the CE's is implemented, indicate that pinning of defects (at w ^ pb ) and generation of 
defects (at w = pbo) by this (phase) disorder in the periodic potential forms the mechanism 
for the smoothening of f c . 

In conclusion, we have shown that the critical current of incommensurate vortex flow 
channels is drastically reduced as compared to the expected value for commensurate chan- 
nels. The dynamic behavior exhibited a crossover from defect motion to coherent flow. We 
described these transport characteristics analytically for the case of a single vortex row. 

We would like to thank A.E. Koshelev, A. van Otterlo and T. Droese for stimulating 
discussions. 
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FIGURES 

FIG. 1. (a) f-v characteristics for vortex channels with w ~ bo an d Ax = 0. Symbols are 
simulated data and drawn lines are obtained with Eq. (|5|). The inset shows an expanded view 
of the small velocity regime. Left upper corner shows the channel geometry with mobile vortices 
displayed by open symbols, (b) Vortex density p(x) = l/{w{xi + \ — Xj)) along the channel for 
w = 0.976o: (o) GS as well as a snapshot of the moving row in the low mobility regime (/ = 0.01), 
(•) snapshot at large drive (/ = 0.1). The drawn line represents p(x) = (w(h(ip + a) — h{ip) + a))~ 1 
at v{f = 0.01). 

FIG. 2. Shear strength f c versus channel width for a velocity criterion of 5 x 10~ 4 and 
Ax = Ax(w). Open symbols are experimental data of artificial vortex channels at the same 
criterion [||. The dashed line represents the continuum result: f c = fibo/w. The inset displays 
typical simulated f-v characteristics. Note the hysteresis for w = 4.54&o and w = 3.76&o- 

FIG. 3. (a) Two snapshots of the channel flow for w = 4.08&o (case I) in the low mobility 
regime (f=0.007, see inset to Fig. 2). For clarity the amplitude of vortex displacements from the 
average y-coordinate of a row has been enlarged by a factor 3.5. Snapshot 2 (o) is taken after an 
average displacement of O.lao with respect to snapshot 1 (•), whereas the defect displacement is 
about 6ao- The rows are labeled with their interstitial density and average velocity. Defects are 
indicated by arrows, (b) (o) defect structure at w/bo = 4.54 (case II), small dots represent the 
trajectories in the plastic regime at / = 0.002. 
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